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Abstract 


The  formali^sm,  developed  by  Kato,  which  gives  upper  and  lower 

bounds  on  the  phase  shift  for  the  scattering  of  a  particle  from  a 

center  of  force   is   shown  to  have  useful  applications  when  the 

scatterer  is  a  compound  system.     In  particular,    the  problem  of 

low  energy  positron  scattering  from,  atomic  hydrogen,  with  zero 

total  orbital  angular  momentum,    is  studied.      It  is   shown  that  at 

zero  energy  the  ordinary  variational  calculation,  which  ignores 

second  order  contributions,  provides  a  bound  on  the  scattering 

length,   from  which  a  bound  on  the  cross  section  is  deduced.     For 

non-zero  energies  a  bound  on  the  phase  shift  may  similarly  be 

obtained,  but  for  a  fictitious  problem  \-ilth  cut-off  potentials. 

If  the  energy  is  sufficiently  small   (less  than  3  e.v.,   say)   the 

error  thus  incurred  is  expected  to  be  negligibly  small.     Numerical 

calculations  performed  at  k  =  0  and  ka     =0.2  lead  to   the  result 

o 

that  the  effects  of  polarization  are  large  enough  to  cause   the 
positjxsn  to  be  on  the  whole  attracted  to  the  hydrogen  atom. 

It  is  shown,   independently  of  the  Kato  formalism,   that   in 
general  a  bound  on  the  scattering  length  may  be  obtained  from  the 
ordinary  variational   calculation  provided  no  composite  bound 
state  exists  for  the  system. 


A  preliminary  report  was   given  at  the  Washington  Meeting  Ox   the 
American  Physical  Society  in  April,   1958    [Bull.   Am.    Phys.  Soc, 
Ser.   II,  3,  171  (1958)] 
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1.     Introduction 

It  is  somewhat  disconcerting  that  more  than  thirty  years  after  the 
advent  of  quantum  theory  there  is  in  common  use  no   general  method 
that  may  be  applied  v/ith  confidence  to   the   calculation  of  the  cross 
section  for  scattering  by  a  compound  system.      This   is  true  even  for 
compound  system.s  consisting  of  only  two   particles.     In  the  few  cases 
in  which  the  calculation  is  on  finn  theoretical  grounds,   such  as   the 
low  energy  scattering  of  neutrons  from  chemically  bound  protons'-  ■*  j 
the  reliability  of  the  calculation  depends  on  some  quite  special 
properties  of  the  forces  and  of  the  scattering  system. 

It  is  of  course  true  that  there  are  a  few  scattering  problems 
involving  compound  scatterers  which  do  not  possess  any  special  properties 
for  which  the   calculation  can  be  accepted  with  a  high  degree  of 
confidence.     An  example  is  the  scattering  of  electrons  by  hydrogen 
atoms '- -^ ' '•-^  ( fo  r  restricted  ranges  of  scattering  energies). 
The  fact  remains,  nervertheless,   that  these  calculations  do  not 
give  any  indication  cf  their  reliabilityj  they  do  not  provide  any 
estimate  of  the  error  nor,   contrarj'  to  the  situation  in  the  evaluation 
of  the  energy  of  the  ground  state,  do  they  provide  an  upper  or  a  lower 
bound  on  the   cross   section.     Rather,   our  confidence  is  based  on  the 
fact  that  variational   calculations  performed  with  a  variety  of  trial 
functions  lead  to  results  which  do   not  differ    appreciably  among 
themselves,   nor  with  a   result  based  on  the  accurately  determined  bound 
state  wave  function  of   H"*-^,  The  use  of  a  variational  method  is  in 
itself  by  no  means  sufficient  to  ensure  any  accuracy  whatever.     Since 
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the  effects  of  polarization  can  be  profound,   it  may  be  quite  difficult 
to  determine  a  suitable  form  for  the  trial  function,  and  of  course  even 
in  a  variational  calculation  a  poor  trial  function  generally  gives 
a  poor  result.     Further,   the  inclusion  of  more  parameters  in  the  trial 
function  need  not  improve  the  result.     In  applications  of  the  Kohn 
and  Hulthen  forms  of  the  variational  principle  certain  consistency 
criteria  are  frequently  used  to  check  the  validity  of  the  calculation'--^. 

It  is  Quite  clear,  however,   that  these  criteria  are  by  no  means 
completely  satisfactory.     The  knowledge   that  the  consistency  ratio   is 
very  different  from  unity  is,   to  be  sure,   a  useful  if  negative  piece 
of  information;   the  difficvJ.ty  is  that  a  consistency  ratio  close  to 
unity  may  iiriply  notl'dng  more  than  an  accurate  calculation  within  the 
limitations,  however  severe  and  inappropriate,  of  the  assumed  form  of 
the  wave  function.     Thus,  Bransden,  Dalgamo,   John,   and  Seaton"-  -"   offer 
the  example  of  a  trial  function  for  the  electron-hydrogen  scattering 
problem  in  the  form  of  a  properly  symmetrized  product  of  one-particle 
functions.     The  assumed  form  gives  rise  to  certain  (approximate) 
one-body  equations.     If  the  one-particle  functions  are  determined 
exactly  as  the  solutions  of  these  equations,  perfect  consistency  will 
be  achieved  even  if  the  resultant  trial  function  bears  no  relationship 
to  the  time  function,     h  still  simpler  example  of  the  failure  of  the 
consistency  criteria  is  found  in  the  present  e  H  problem.     A  trial 
function  involving  no   polarization  of  the  hydrogen  atom  will  lead  to 
perfect  consistency  if  the  positron  wave  function  is  the  exact  solution 
of  the  static  problem,   and  yet  polarization  plays  a  crucial  role.     The 
weakness  of  the  consistency  checks  is  well  recognized,  but  though  all 
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agree  with  s   previous  author"-  -'  that  "a  foolish  consistency  is  the 
hobgoblin  of  little  minds",  these  checks  have  nevertheless  been 
frequently  applied  simply  because  no  better  standards  have  been  available. 

In  the  case  of  scattering  by  a  static  central  potential,  where  the 
cross  section,  cr,  consists  of  a  sum  of  partial  cross  sections,  <X, 
each  of  which  depends  upon  just  ony  real  number,  the  phase  shift Tn,, 
a  number  of  methods  have  been  devised  for  determining  upper  and  lower 
bounds  on  cot  "Vi^  and  hence  on  cTt   •  ^y  ^^^  *^^  '^^'^  general  of 

re"! 

these  methods  is  the  one  due  to  Kato^-'.  The  method  was  applied  by 

Kato  to  the  case  of  L  =  0  with  rather  striking  results.  It  has  also 

[9! 

been  applied  to  higher  angular  momenta '■-^  J  the  calculations  are  then 

more  cumbersome,  but  useful  results  can  still  be  obtained.  Some 
improvements  in  the  formalism  itself  have  also  been  effected,  including 
the  minimization  of  the  amount  of  a  priori  knowledge  of  "^  which  is 

Fiol 

required  before  the  method  can  be  applied '-  -"  j  this  will  generally  be 
essential  for  the  compound  scatterer  where  (except  for  zero  energy)  one 
ordinarily  has  no  a  priori  knowledge  of  "h-  whatever. 

It  is  to  be  noted  that  the  Kato  formalism  is  by  no  means  restricted 
to  quantum,  mechanical  problems;  it  is  applicable  to  a  wide  variety 
of  problems,  including,  for  example,  the  study  of  electromagnetic 
waveguides '-  -^ . 

It  is  to  be  expected  that  some  difficulties  will  arise  in  an  attempt 
to  generalize  the  Kato  method  to  scattering  by  a  compound  system.  These 
difficulties  will  be  discussed  in  Section  2  where  it  will  also  be  shoxm 
that  there  exist  interesting  circumstances  under  which  they  can  be 
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circumvented.  It  is  further  to  be  anticipated  that  a  calculation  which 
furnishes  one  or  both  bounds  will  be  more  complicated  than  a  variational 
calculation.  While  this  is  generally  the  case,  there  are  situations 
in  which  an  ordinary  variational  calculation  gives  a  bound. 

Of  the  problems  to  which  the  Kato  method  can  be  applied  we  have 
thought  it  worthwhile  to  consider  first  the  scattering  of  positrons 

ri2i  ri3i 

by  hydrogen  atoms'-  -^ '  ^    -^  .   This  problem  is  of  some  interest  itself  j 
also,  the  methods  involved  can  be  extended  in  a  natural  way  to  the  study 
of  e  on  He  where  a  direct  check  between  theory  and  experiment  is 
possible^"'"^-'*^  -'''--'.  The  problem  has  the  additional  virtue  of 
providing  an  opportunity  for  studying  the  effect  of  "  specific 
polarization"  ,  i.e.,  the  distortion  of  the  scattering  system  which  is 
present  in  the  wave  function  whether  or  not  its  antisymmetrization  is 
required.  Application  of  the  Pauli  principle  automatically  introduces 
some  distortion,  called  "  exchange  polarization"  ,  and  in  the  case  of 
electron  scattering  hy  hydrogen,  for  example,  it  is  rather  difficult 
to  separate  out  these  two  "  polarization"  effects^  -' .  However,  the 
primary  reason  for  considering  this  problem  first  is  its  simplicity. 
To  begin  with,  it  is  an  atcmdc  problem  which,  being  free  from  the 
additional  complexities  of  tensor  forces  and  space  exchange,  presents 
a  more  convenient  starting  point  for  this  type  of  investigation.  The 
absence  of  the  Pauli  principle  and  the  fact  that  no  three  body  bound 
state  exists  further  simplify  the  analysis.  The  "  infinite"  mass  of 
the  proton  as  compared  to  the  electron  or  positron  mass  is  also  a 
simplification  since  the  complexities  that  sometimes  arise  from  the 
separation  of  the  center  of  mass  cannot  occur.  Of  course,  the  possibility 
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of  positron-electron  annihilation  exists.     If,  however,  we  assume  that 
the  annihilation  cross  section,  cr^is  given  byi-    -I 

(l.l)  cr    =  (very  small  constant)  v"  , 

a 

where  v  is  the  speed  of  the  positron,   it  follows  that  only  at  extremely 
small  energies  will  this  process  compete  with  elastic  scattering. 
(Actually,   the  above  formula  ignores   the  velocity  distribution  of 
the  electron  in  the  atom,   as  well  as  the  electrostatic  interaction 
between  positron  and  electron.     These  effects  might  well  reduce  the 
annihilation  cross  section.)     Our  subsequent   "  zero"  energy  results 
may  therefore  be  interpreted  as  extrapolations  from  some  small 
non-zero  energy.     The  situation  is  quite  similar  to   that  which  obtains 
in  the  case  of   "  zero"    energy  neutron-proton  scattering,  where  the 
competing  capture  process  obeys  a  v~     law  and  is  significant  only  for 
the  lowest  energies.     The  only  unpleasant  feature  associated  with  the 
e*H  problem  is  the  presence  of  the  long  range  Coulomb  potential. 
The  numerical  calculations  are  discussed  in  Section  U. 

Some  general  remarks  should  be  made  about  the  utility  of  theoretical 
calculations  of  the  type  under  consideration.     We  note,   first,   that 
to  our  knowledge  no  machine  calculations  have  been  performed  which 
rigorously  treat  the  scattering  by  a  compound  system.     Further,   in 
electron  scattering  by  compound  systems,   the  region  of  greatest  interest 
is  the  low  energy  region  which  is  at  the  present  time  relatively 
inaccessible  to  experimental  investigation.     It  is  true  that  one  is 
primarily  interested  in  electron  scattering  by  oxygen  and  nitrogen, 
due  to  its  importance  in  ionospheric  processes,   and  the  present  methods 
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cannot  yet  be  applied  to  such  a  complicated  case.  On  the  other  hand, 

a  theoretical  analysis  of  low  energy  sin-^let  and  triplet  scattering 

of  electrons  by  hydrogen  atoms  is  possible,  (work  on  this  is  in  progress), 

and  a  close  estimation  of  the  cross  section  for  this  process  should 

be  useful  in  the  determination  of  the  accuracy  of  various  approximations; 

those  v;hich  are  found  to  be  reasonably  accurate  can  then  be  made  with 

a  considerably  greater  degree  of  confidence  in  the  more  complicated 

problems.  Finally,  in  the  low  energy  nuclear  scattering  problems 

(neutrons  on  deuterons,  for  example),  where  the  experiments  can  and 

have  been  performed,  a  theoretical  analysis  is  still  required  for 

in  these  problems  one  is  primarily  interested  in  deducing  the  potential 

from  the  scattering  data, 

2.  The  Kato  Method 

A.  The  True  Problem 

To  facilitate  subsequent  discussion  we  define  a  problem  with 
somewhat  specialized  conditions.  We  consider  the  scattering  of  a 
particle  (particle  number  l)  of  mass  mj  from  a  system  in  its  ground 
state  consisting  of  two  particles  of  masses  M  and  m-  respectively. 
Vie   take  m..  =  nu  =  m  and  choose  M  »  ni  so  that  the  scattering  system  is 
effectively  a  particle  (particle  number  2)  of  mass  m  in  a  field  of  force. 
Particles  1  and  2  are  assumed  to  be  distinguishable  so  that  we  need  not 
concern  ourselves  with  the  Pauli  principle.  Let  E^  denote  the 
ground  state  energy  of  the  second  particle  in  the  field  of  force; 
AEq  is  its  minimum,  excitation  energy.  Correspondingly,  if  the  first 
particle  can  be  bound  to  the  center  of  force,  let  E-,  denote  the  ground 
state  energy.  If  particles  1  and  2  can  by  themselves  form  one  or  more 
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bound  states  we  denote  their  .-ground  state  energy  by   E   .  Finally, 
we  let  T-|  Isbel  the  initial  kinetic  enercy  of  the  incident  particle. 
The  discussion  is  limited  to  scattering  states  of  zero  total  orbital 
angular  mcrrentuin. 

The  Kato  method  is  applicable  only  if  the  scattering  problem  is 
such  that  the  various  angular  moments  are  not  coupled  and  if  the 
scattering  for  a  given  angular  momentum  can  be  completely  characterized 
by  one  real  phase  shift.  For  the  scattering  of  a  particle  by  a 
center  of  force,  both  these  conditions  are  realized  if  the  potential 
is  spherically  sjiranetric.  In  the  many  body  scattering  problem  we  must 
not  only  have  potentials  which  are  independent  of  the  spatial 
orientation  of  the  system  as  a  whole,  to  allov;  a  partial  wave  analysis, 
but  also  T  must  be  restricted  in  magnitude  if  the  characterization  by 
one  real  phase  shift  is  to  be  possible.  In  particular. 


(1)  Excitation  or  ionization  must  not  be  possible;  i.e.,  we  must 

have  T^  <  A  ^2  * 

(2)  Exchange  must  not  be  possible,   requiring    ISp    I   >   '^le''   ^^'^ 
further  T^  <   \E^^\    -    \e^J. 

(3)  Pick-up  must  not  be  possible,    requiring    Is^^l   >   l%2<^^    ^"^ 

^1-=^   l^2gl   -    1^12gl- 

In  the  case  of  positron  scattering  from  hydrogen  atoms,  we  have 

\%\  =i-'A   ^AE^-^lEgJ    ,      V^  =  i^2g   • 

Exchange  is  not  possible  and  conditions  (l)  and  (3)  above  require 

that  T,  be  less  than  -pJE   |  and  -^JE-  |  ,  respectively.  It  is  necessary 
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therefore,  that  the  initial  positron  energy  be  less  than  6.8  e.v. 
in  order  that  one  real  phase  shift  completely  characterize  the  L  =  C 
scattering. 

We  return  now  to  the  more  general  problem  and  assume  that  conditions 
(l),  (2)  and  (3)  are  satisfied.  The  differential  equation  to  be  solved 
is 

(2.1)  (2in/fe^)(E  -  H)  \ir^,T^tr^^)/r^   -  (l/r^)C«^(r^,r2,rj^2^  '  °   * 

where  H  is  the  Hamiltonian  of  the  system  and  E  =  E-  +  T-  is  the  total 
energy.  Here  r^  =  |r^|,  Pg  =  ^2!,  and  r^2  "^  l^j  '  ^2^  >   where  r^ 
and  r-  are  the  position  vectors  cf  the  two  particles  measured  from  the 
center  of  force.  Defining  e  =  (2mAi^)|E2  I,  k^  =  (2m^  )T^,  and 
W  =  -(2m/fi  )V  for  all  three  "  potentials"  W,,  W^  and  W^g*  ^^   ^^^® 

(•?    n\         P  -        ^      ^     JL.        3    ^  2  _B_  ^  ,J^   .   _1_  ^   e  ,,     2^  3 

^'•'^  ^"i;r?    ;7  ^ '2  ^^^    r2^^^^'-^^~p 

.W^(r^).W2(r2)-W^2W2^  -  k^  -  e^  , 

where  p  is  the  cosine  of  the  angle  between  r-,   and  r2.     It  is  sometimes 
more   convenient  to  write  the  total  kinetic  energy  operator,  T,   as 

(2.3)  .j„A>^i=  ^  J-  ^i'^*T^  -/  J: 

r^  1  1       rg         ^  -i 

2     ar,„     ^12     ar,o  r,r,o  bj^hr 


r^^-         12       ^^     °^12  '■L'12  "^n-^12 

(r  2+r     ^-r  M  2 


^2^12  ^^2^^12 
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The  exect  solution,  iL,  must  satisfj"-  the  boundary  conditions 

^  _>  Kr^)  cos(kr^+e)  +  cotC'vr-e)  3in(kr^+&)  ,  as  Tt^ ->oo 

(2.U)  Uq  ->  0  ,      as  r^ 

Uq  =  0      at  r^  =  0 

where  R(r„)  is  the  exact  ground  state  wave  function  for  the  bound 
particle,  normalized  such  that 


^2  —'  00 


(2.5)     2  J     r^^dr^  R^(r2)  =  1. 


The  normalization  parameter,  ©,   satisfies  0  <  &  <  n  but  is  otherwise 
arbitrary.     Tlie  true  phase  shift, "fT  ,   is  determined  by  (2.1),    (2.2) 
and  (2, It).     (We  have  here,   and  in  the  following,  dropped  the  subscript 
zero  on'vi  ,) 

We  now  introduce  a  trial  function,  UgCrn*  r^,  ^■]2^ >  which  satisfies 
the  same  boundary  conditions  as  iL  but  with  the  true  phase  shift,  7)   , 
replaced  by  a  trial  phase  shift,  >)  .     Note  that  this  still  involves 
the  exact  wave  function,  R{rJ).     The  identity 

(2.6)  k  cot(>i-9)   =  k  cot(v|-e)-      u^Xugdr  +       w^i^w^dT 

2 

then  carries  over  from  the  one  body  case  where  now  df  =  dr,   r^     drp  dp 

and  where 


(2.7)     ^Q^^2.f^2'^12^  "■  '%^^1'^2*^12^"^^^1»^2'^12^* 
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Since  w^j  is  a  first  order  term,   Eq.    (2.6)   constitutes  a  varistional 
principle  for  cot   (y)-9)  upon  dropping  the  last  term. 

The  requirement  that  the  asymptotic  form  of  u„  be  proportional 
to  the  exact  wave  function  R(r_)   is  evidently  extremely  restrictive. 
It  is,  however,  necessary  in  order  that  the  decomposition 

(2.8)  ^^\  d*^    =    Xu^oCugdr-J  Wg^w^dt 

may  be  used.  If  the  requirement  is  dropped  the  integrals  on  the  right 
hand  side  of  the  above  equation  will  not  separately  exist,  each  of 
them  diverging  for  large  r,.  (The  sum  of  the  two  integrals  is  of  course 
still  finite.)  The  required  knowledge  of  the  ground  state  wave 
function  formally  eliminates  the  possibility  of  applying  the  Kato 
method  to  practicalDy  all  scattering  systems  consisting  of  more  than 
two  particles.  This  difficulty,  present  in  the  ordinary  variational 
problem  as  well,  can  fortunately  be  partially  circumvented.  Assuming 
a  knowledge  of  the  ground  state  wave  function  we  introduce  into  the 
appropriate  integral  expressions  a  trial  function  which  is  the  sum 
of  an  "inside  wave  function"  ^  '-J  arid  a  term  which  provides  the  correct 
asymptotic  fonaj  the  operation  of  X  on  the  asjTnptotic  form  is  carried 
out  formally.  At  this  stage,  we  can  replace  the  ground  state  wave 
function  by  some  approximation  to  it  without  introducing  any 
infinities.  The  error  involved  will  be  reproduced  in  our  estimate  of 
(or  bound  on)  k  cot  (y[-©)  to  first  order,  i.e.,  this  error  is  not 
reduced  by  the  fact  that  a  variational  formulation  is  being  utilized. 
Further,  the  question  as  to  whether  the  bound  had  been  preserved  would 
have  to  be  examined.  However,  for  the  cases  for  which  the  ground  state 
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wave  function  is  knoi'm  with  considerable  precision,  (He  for  example), 
the  replacement  should  be  irrelevant. 

B.  The  Associated  Eigenvalue  Problem. 
Consider  the  equation 

(2.9)    (£0^  +  i^g  p  0Q  =  0 

where  p  is  a  non-negative  function  to  be  chosen  such  that  the  scattering 
problem  determined  by  Eq.    (2.9)  and  by  appropriate  boundary  conditions 
is  characterized  for  each  value  of  ii,_  by  one  real  phase  shift,  6(iJ.g). 


Then  there  exists  an  infinite  set  of  discrete  eigenvalues  u    a, 

corresponding  eigenfunctions  0     r^  normalized  by 

n," 


ind 


(2.10)         k"^ 


y 


0„  a  0„  a  P  dr  =     6         ,  in,n  =  0,  ±1,  ±2,    ... 

n,w     m,w  nm 


where  5       is   the  Kronecker  5-symbol,   such  that 
ran 

0    ^  ->  const.  R(rr,)  sin(kr,  +  &  +  nn),     as  r,   ->  oo 
n,fcy  ^1  -L 

(2.11)  0^  ^  ->  0  ,  as  Tg  ->  00 
^n,e  =  °           '  for  r^  =  0   . 

We  denote  the  smallest  positive  eigenvalue  of  the  set  p.  „  ^y  ^■o  ^""^  ^^'^ 
smallest  (in  magnitude)  negative  eigenvalue  by  -Pq.  The  Kato  inequalities 
are 

(2.12)  ^^'^  j  (JCuq)^  p"-^dr<  I  WqXwq  d-t  <  Pq"^  J  OCu^)^  p'-'-dr    . 

This  formal  development  follows  that  for  the  one  body  problem, 
mutatis  mutandis >  (Some  questions  of  completeness  are  discussed  in 
Appendix  B).   However,  the  actual  restrictions  on  p,  referred  to  above. 
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involve  considerations  peculiar  to  the  many  body  problem.  We  observe  first 
a  restriction  which  does  carry  over  from  the  one  body  case,  namely,  that 
p  must  vanish  faster  than  l/r-,  as  r,  becomes  infinite  in  order  that 
normalization  according  to  Eq.  (2,10)  be  possible  and  the  phase  shift, 
6(|j.),  be  defined.  The  choice  p  =  p(rp),  for  example,  is  eliminated'"''. 
To  eliminate  pick-up  in  the  associated  eigenvalue  problem,  which  is 
defined  for  an  infinitely  broad  range  of  potential  strengths,  there  are 
two  alternatives. 

a)  Choose  p  to  be  independent  of  r, „,  at  least  for  r,  and/or  r^  ->  oo. 
Then  E^.  will  be  independent  of  \<.. 

b)  Imagine  the  bound  particle  to  be  enclosed  in  a  box  with  perfectly 
rigid  walls.  The  size  of  the  box,  although  finite,  may  be  large  enough 
so  that  the  deviation  from  the  real  problem  is  negligibly  small. 

The  exchange  process  may  occur  in  the  associated  problem  even  if 
it  is  absent  in  the  real  problem.  Again  there  are  two  alternative 
solutions. 

a"*)  Choose  p  to  be  independent  of  r, ,  at  least  for  Tp  ->  oo.  Then 
E_  will  be  independent  of  |i. 

b')  Put  the  bound  particle  in  a  box,  as  in  b)  above. 
If  this  latter  suggestion  is  adopted  we  require 

(2.13)   r^p(r^,r2,r^2^  ->  0  ,    as  r^  ->  oo, 
but  p  is  otherwise  arbitrary. 

■it- 

The  theorem  which  states  that  the  phase  shift  for  scattering  by  a  center 

of  force  is  a  monotonic  function  of  the  potential  strength  for  a  potential 

of  constant  sign  carries  over  to  the  many  body  case  provided  that  the  unperturbed 

energy  of  the  compound  s3''stem  remains  fixed  as  the  potential  strength  varies. 

Thus,  another  reason  for  demanding  that  p  vanish  with  large  r-,  is  that  the 

required  monotonicity  of  6{\x)   with  jj.  will  not  be  preserved  if  p  becomes 

some  function  of  r2  as  r-,  ->  oo.  Furthermore,  for  large  negative  values  of 
H,  particle  2  would  no  longer  be  bound.  The  boundary  conditions,  Eq.  (2.11), 

would  then  be  violated. 
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3.  Application  to  Positron  Scattering  From  Hydrogen. 

Kato  has  shown  for  the  static  central  potential  that  if  the  potential 
cuts  off  outside  a  region  of  radius  a  then  p^  =  oo  if,  in  addition, 
the  following  conditions  are  satisfied 


(1)  ka  <  TT-© 

(2)  v[  <&<n 


^  /,   9   -1 

(ivu)  p  dt   .  The  analogous  result  holds  in  the 


This  is  a  particularly  interesting  case  since  then  a  lower  bound  may  be 
obtained  on  the  phase  shift  without  the  necessity  of  evaluating  the  more 
difficult  integral, 
three  body  problem  considered  in  Section  2.  In  this  case  the  required 
proof  that  5(-oo)  =  -ka  (where  p  =  0  for  r^  >  a)  is  slightly  more  involved. 
Vfe  may  write  the  solution  for  r^  >  a  as 

(3.1)  0   =  const.  R(r2)  sin(kr^  +  6(ij.)) 

n,L  J  * 

where  R  T(r„)  are  the  excited  inner  particle  bound  states  and  the 
n,L  <£ 

continuum  solutions,  hJ   are  the  spherical  Hankel  functions  and  P  are 
the  Legendre  polynomials j  k  is  defined  by 

(3.2)  k^-Cg^k^.    e^     . 


Fo 


r  10.  =  -00  we  must  have  ^(r^   =  a)   =0.     From  the  orthogonality  of  R  with 


the  functions  R     -  we  have 
n,L 


(3.3)       sin(ka  +  6(-oo))  =  0       , 
from  which  the  result  follows. 
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The  many  body  problem  is  generally  complicated  by  the  difficulty 
that  one  cannot   get  even  crude  bounds  onvj   ,   so   that  condition  (2)   cannot 
be  proved,    (assuming  that   it  is  true).     We  have  recourse  here  to   the  use 
of  a    "  conditional  inequality"  *-    -• .     We  simply  assume  that  condition  (2) 
is  satisfied;    the  subsequent    "  conditional  inequality"    for  k  cot(>}  -  9  ) 
may  or  may  not  be  valid.     However,  the  bound  onT?  deduced  from  the 
"  conditional  inequality"    on  k  cot(ri"-&)  wil]   be  valid  whether  or  not 
the  bound  on  k  cot("»7  -&)   is  valid,   and  therefore  whether  or  not  condition 
(2)   is  in  fact  satisfied. 

We  therefore  consider  the  following  set  of  potentials 

(3.1i)         W2(r2)   =  SCa^r^)"^         ,  0  <  r^  <  b 

=   c  ,  b<r2<b  +  d 

=  -(X>  f  r_>b  +  d 

(3.5)         W^(r^)  .W^2(^12^  =  -^^  Vl^"^     ^2(a^r3^2^-l  ,     r^  ^  a 

=     0  r^  >  a 

Here  a  is  the  Bohr  radius.  As  a  and  b  become  infinite  the  above 
o  —    — 

potentials  approach  those  of  the  positron- hydro gen  system.  As  indicated 
above,  b  may  be  chosen  large  enough  so  that  the  presence  of  the  infinitely 
repulsive  potential  barrier  has  no  effect  on  the  real  problem",  although 
it  plays  an  important  role  in  lim.iting  the  number  of  allowable  processes 
which  can  take  place  in  the  associated  problem.  Furthennore,  for  k 

Strictly  speaking,  one  should  study  the  integral,   rudou^dTi  ,  in  the 

limit  of  large  b  ,  with  the  parameters  c  and  d  always  adjusted  to  give 

the  correct  hydrogenic  energy,  e  j  the  error  involved  in  changing  the 

S 
potential  should  then  be  shown  to  be  arbitrarily  small  for  b  sufficiently 

large . 
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sufficiently small,  a  may  be  chosen  large  enough  so  that  the  neglect 
of  the  potential  which  may  exist  beyond  r-i  =  a  introduces  a  negligible 
error  in  the  phase  shift.  If  k  is  zero,  we  may  let  a  become  infinite, 
in  which  case  we  introduce  no  distortion  of  the  positron-hyxiiK?gen 
potentials.  In  this  case,  the  validity  of  the  second  condition  stated 
above,  namely ">?  <  6  <  n,  is  deduced  by  assuming  that  no  three-body 
bound  state  exists  for  the  positron-hydrogen  system  (none  has  been 
found  either  experimentally  or  theoretically) '-  -' .  We  then  apply  a 
theorem  which  is  the  generalization  to  many  body  systems  of  a  theorem  proved 
by  Levinson^  '-'  for  the  scattering  from  a  static  potential.  It  states 
that  when  the  exclusion  principle  is  not  in  effect  the  phase  shift 
for  zero  energy  scattering  is  nn  where  n  is  the  number  of  bound  states 
of  the  composite  system'-  -^ .  From  this  we  conclude  that  17  =  0  for  k  =  0. 
A  rigorous  proof  of  the  above  theorem  has  not  been  given  even  for 
potentials  which  fall  off  rapidly.  A   proof  for  the  slowly  decaying 
Coulomb  forces  will  be  even  more  difficult.  We  nevertheless  apply  the 
theorem  with  some  confidence,  for  its  generalization, which  applies  to 

electron  scattering  from  atoms,  has  proved  to  be  correct  in  all  cases 

fPl! 
where  it  has  been  checked"-  -■ .  That  the  theorem  should  be  applicable  to 

the  Coulomb  case  (assuming  it  applies  to  rapidly  decaying  potentials) 

is  perhaps  made  more  reasonable  by  the  observation  that  the  "  effective 

interaction"  of  the  positron  or  electron  with  the  atomic  system  falls 

off  fairly  rapidly.  For  e"*"  on  H,  e.g.,  it  goes  as  Vr  ,  Actually,  as 

we  will  show  in  Section  ^,   the  results  we  have  obtained  do  not  depend  on 

the  validity  of  the  above  mentioned  theorem. 

With  the  choice  &  =   n/2  and  -A  s  lim  "yj/k   ,  -A  =  lim  l^/k, 

k->o   ••  k^o  ^ 

u  «  lijB    *'(i/2)n^»  *®  ^^"^^  **^  following  bound  on  the  scattering  length,  A, 
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UCM  u     d^ 


(3.6)  A     <    A  - 

J 

where  u  has  the   asymptotic  form 


One  of  the  strong  features  of  the  Kato  method  is  that  in  general 
it  provides  a  different  bound  for  each  choice  of  9.     This  is  not  true 
at  k  =  C,  however,  where  the   identical  result  is  obtained  for  any  & 
other  than  9=0,     The  choice  &  =  0  will  lead  to  a  valid  bound  on  A 
only  if  "Yi   approaches  zero  from  below  as  k  goes  to   zero,  i.e,,   only  if 
A   is  positive.      In  this  case  we  have 

(3.8)         A   -^  >     A"^+    A"^       tlcTu  dr 


A   -1  >     A-^  .    A-2  J  {IjTu 


If  u  is  sufficiently  close  to  the  true  function  such  that  A  and  A 
are  of  the  same  sign  and  such  that    |ll«iA|,  where 

i^jsC^dtf 

then  Eq.   (3.8)  may  be  rewritten  as 

(3.8)      A<  A  -  I  +  A~^  i^  +  ...,  |i|«|aI. 

It  then  follows  (since  A  >  O)  that  Eq.  (3.6)  is  superior  to  Eq.  (3.8), 
Therefore,  for  the  special  case  of  zero  energy  scattering,  where  no 
composite  bound  state  exists,  the  problem  of  the  optimum  choice  of  © 
is  particularly  simple.  One  should  use  Eq.  (3.6),  which  follows  from 
any  choice  of  9  satisfying  0  <  &  <  n  , 
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k.     Numerical  Calculation 

A.  No  Polarization  Approximation,  k  =  0. 

To  begin  with,  we  consider  the  class  of  trial  functions  of  the 
form 

(it.l)   u(r^,r2,r^2^  ^   g(r2)f(r^)  . 

The  required  asymptotic  form  fixes   g(r,)   as 


gv^2 


We  also  have 


(rj   =R(rJ   =   (2/a  3)1/2  e-^2Ao. 


(U.2)  f(o)   =  0 

f(r, )  ->  A-  r-     ,  as  r^  ->  oo. 

A   trial  function  in  this  product  form  corresponds  to  the  approximation 
(variously  referred  to  as  the  no  polarization  approximation,  the  static 
approximation,   and  the  one  body  approximation)  in  which  the  hydrogen 
atom  is  not  polarized  by  the  incoming    positron.  We  substitute   the  above 
expression  for  u  into  Eq.    (3.6),  and  integrate  over  r^  and  .p.     Dropping 
the  subscript  on  r^   we  find 

fOO 


(Li.3)       A     <    A  -         f(r)  <C    f(r)di 


where 

(b.U)  J:     =     -^     +  W  (r) 

dr 

and 

-2r/a 
(U.5)  W  (r)  =  -(2/a^)  e  °   (l/a^  +  l/r) 
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so  that  in  the  static  approximation  the  positron  is  repelled  by  the 
hydrogen  atom.      The  best  possible   choice  for  f(r)  would  be  T  (r),   the 
exact  wave   function  for  scattering  by  a  fixed  potential,  W  (r).     This 
T-jill  now  be   shown  to   follow  from  an  application  of  the  Kato  method  to 
this  one  body  problem'.     Since  no  bound  state  exists  for  the  potential 


W  (r)   +  np(r)  with  [i  negative  we  conclude   that 


s 


(I1.6)  6(ti)   =0  ^i  <  0 


and  therefore  tha 

we  have 

,00 


t  Pa /A^  =    00.     Thus,   for  any  trial   function,   f(r). 


(U.7)         A     <    A    - 


y 


f(r)cC     f(r)d3 


where  T.     is  the  scattering?  length  for  the  static  problem.  Since  the 

equality  holds  only  when  f(r)  =  f  (r),  the  statement  above  has  been 

verified.  From  Eq.  (li.3)  with  f(r)  =  f  (r)  we  find  that 

s 


(U.8)        I<I 
—     s 


In  point  of  fact  one   cannot  solve  for  f  (r).     However,   since  we 

are  dealing  here  with  the  scattering  by  a  static  potential,   at  zero 

angular  momentum  and  zero  energj'',   it  is  a  trivial  matter  to   find  an  upper 

bound  on  A     and  hence  on  A.     V7e  will  also  find  a  lower  bound  on  A     in 
3  s 

order  to  measure   the  accuracy  of  our  determination  of  k     and  to  exhibit 


It  should  be  emphasized  that  in  the  remainder  of  sub-section  A  various 
quantities,  and  in  particular,   pv.  and  a^,   refer  to  this   (static) one  body 
problem  and  not  to  the  tine  problem. 
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once  again  the  power  of  the  Kato  method  for  this  type  of  problem. 

We  have  already  shown  that  P/-,  /A^  =  oo.     Choosing  p(r)  =  -W^(r) 
we  observe  that  at  ,^   is  the  smallest  positive  value  of  iJ.  for  which 
the  "  potential"  W  (r)(l-|i.)  can  support  a  bound  state,  for  the 
appearance  of  a  bound  state  indicates  that  the  phase  shift  has  jumped 
from  zero  to  n.   (T^ctually,  ^i%A^^   =(l/^.)  Since  Wg(r)(l-n,)  isn't 

even  attractive  until  n  =  1,  we  immediately  obtain  °-h/Aj.  >  !•  This 

[22I 
lower  bound  can  be  improved  quite  easily.  Bargmann*-  J  has  shown 

that  a  necessary  condition  for  the  existence  of  a  bound  state  for 

a  potential  U(r)  is 


rOO 

(U.9) 


,00 

I  r  Iu(r)l  dr>  1. 


With  U(r)  =  -(1-|J,)W  (r),  this  leads  to  oa/A^  >  5/3.' 


If  we  choose 


-nr/a         -qr/a 
(U.IO)    f(r)=A(l-e     °)  -r(l+Be     °) 


we  find,  with  n  =  1.5  and  q  =  2,  that 


(U.ll)         0.5762a^<A      <    0.5823a 


where  A  =  0,590[ia  ,  B  =  -0.1819  for  the  upper  bound  calculation  and 

A  =  0.58l6a  ,  B  =  -0.0698  for  the  lower  bound.      (Even  with  B  =  0   the 
o 

bounds  are  within  5  %  of  each  other.)     /*t  the  same  time  we  have 
(U.12)  A  <    0.58233^       . 

A  comparison  of  (l-p.  )W  (r)  with  the  Hulthen  potential  L  -'J   leads  to 
aL  ,  \    >  2.     Although  we  shall  use  this  fact  later  in  a  slightly  different 
connection,  we  prefer  to  base  the  present  result  on  an  application  of  the 
Bargmann  theorem  rather  than  on  this  rare  occurrence  of  an  excellent 
comparison  potential  for  which  solutions  may  be  obtained  analytically. 
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_2 
The  cross  section  at  zero  enersy,   (T(k  =  o),  is  UnA  .  The  above 

result  without  an  accompanying  lower  bound  does  not  exclude  a  large 

negative  value  of  A  and  therefore  does  not  provide  an  upper  or 

lower  bound  on  cCk  =  o). 

It  is  interesting  to  note  that  a  lower  bound  on  A  could  have 

been  determined  without  performing  a  variational  calculation,   (Since 

this  is  a  lower  bound  on  A  ,  it  does  not  provide  any  bound  on  A.) 

s 

It  has  been  shown'-    -^    that  for  two    "potentials"     +W(r)  and  -VJ(r)  with 
phase  shifts 'h"_^  and  "^  ,  i^spectively,  where  W(r)  >  0,  yj^  <  v.  and 
"h     >  -n,   the  following  inequality  holds: 

()ul3)         k  cotv   -  k  cot>)     <     2  k  cof 


,11^  -  k  coth     <     2  k  cotV- 


Here  "rt  „  is  the  value  of  ]n     calculated  in  the  Bom  approximation.     For 
W(r)  we  choose 

-2r/a 
(U.lii)         W(r)  =  -Wg(r)   =   (2/a^)  e  °   (l/a^ -H/r)  , 

the  stctic   "potential"  for  electron  scattering  by  hydrogen.     Since  -V7(r) 
obviously  cannot  support  a  bound  state,   and  since  it  can  be  shown  that 
W(r)   cannot  support  a  bound  state  either,  we  have  ))^     =  )7       =  Oj   for  k  =  0 
Bounds  on  the  L  =  0  phase  shift,  yj     ,   for  the  above  potential  have  been 


Here  we  do  take  rer.ourse  in  a  comparison  with  the  Hulthen  potential. 
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determined.     In  the  limit  of  zero  energy  we  have'-   -^ 
(U.15)  0.10560  a^"-^  <     (1^  cot-n^g)j^^^  <    0.10592  a^"-^ 


Now,   as  k  ->  o,   the  Born  approximation  yields 


1... '  -  v| 


,-00      ^ 

(['..16)  k  cotV)^  „  ->  -  1/  (       r     W  (r)dr  =  l/a, 

I  J-o  wig  (I 


It  therefore  follows  that 


(li.lY)  (k  cofy"_^)j^^^  >  -Sa^"-"-  +  0.10560  a^""'" 


and 


(U.18)  Ig  =  -(k  coth_g);;i^     >    0.5273  a^ 


While  this  result  is  not  as  good  as  the  one  given  previously,   it  was 
obtained  without  constructing  a  trial  function. 

Summarizing  the  results  of  thj.s  sub-section,  we  find  that  an 
analysis  of  the  static  problem  loads  to  an  upper  bound  or  the  tnte 
scattering  length,  A,  but  does  not  provide  any  information  on  a  lower 
bound  for  ~.     Further,  the  least  upper  bound  obtainable  in  the  static 
approximation,   namely  that   given  by  Eq.    (h.3)>   is  positive  so  that  we 
do  not  obtain  an  upper  or  a  lower  bound  on  o'(k=o).     Actually,  as  we 
will  show  in  the  next  sub-section,  T  is  negative  so   that  the  static 
approximation  is  totally  inadequate  in  this  problem. 

B.  Polarization  Considered,  k  =  o. 

A  trial  function  of  the  form  u  =  g(r-,  jr^)  may  take  into  account, 
to  some  extent,  the  polarization  of  the  hydrogen  atom.  As  we  can  see. 
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using  arguments  quite  similar  to  those  amployed  in  sub-section  UA, 
this  type  of  trial  function  leads  to  a  positive  scattering  length. 
Thus,   substituting  g(r, ^r^)   into  Eq.    (3.6)  and  integrating  over  p  we 


find 


;00  ,00 


where 


,2 


(U.20)     £,   =  -^     ^    -K    ^  r/  ^  .    ^(A  .  J.  .  i)  .  e 
*       ar^2  r/     ^^2     2     °^2         ^o  ^2       ^>       ^1  ^ 


Here  r^  is  the  larger  of  r,   and  r^.     Since  no  bound  state  exists  for  the 

operator  k-  ,   none  exists  for  Ju.  .     This  is  a  consequence  of  the  fact  that 

the  ground  state  energy  satisfies  a  minimum  principle  and  the  space 

on  which  cC,    is  defined  is  a   sub-space  of  that  forju.    (The  same  result 

also  follows  directly  from  the  fact  that  l/r     -  l/r,    is   always  negative l-    ^-',) 

Therefore,   the  exact  scattering  length,  TT  ,   for<Jj.    is  bounded  by  the 

t      t 

Kato  inequality  with  Pa  /A  =  oo.  It  then  follows  that  A  represents 
the  best  (i.e.,  minimum)  result  obtainable  with  a  trial  function  which 
is  independent  of  r-,^*  Since  the  "positron"  is  always  repelled  (i.e., 
l/r^  -  l/r^  ^  ^^  ''^^^  nonotonicity  theorem^  shows  that  A  is  itself 
positive,  which  verifies  our  original  statement. 

It  is  clear  that  a  realistic  trial  function  should  depend  on  r-, ,  Tp 


See  footnote  on  page  12. 
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and  r, p.  This  dependence  was  chosen  under  the  physically  plausible 
assumption  that  the  polarization  is  largely  accounted  for  by  the 
virtual  formation  of  positronium.  Since  the  latter  has  a  ground  state 
function  proportional  to  exp(-l/2  T^-y^hS)   the  tenri  C  expC-tr-jp/^Q  -  ^^^2/3^)  - 
the  additional  factor  being  reauired  for  convergence  -  was  added  to  a  function 
similar  in  form  to  the  trial  function,  Eq.  (Ii.lO),  which  was  so 
successful  in  the  static  calculation.  VJe  have  chosen 


(U.21)  '^{r^,r^,T^^   =  (2/a^3)l/2L "'2  ^o  ^_ 


X       2'  o 
)  -  r^e 


-qr/a^  -  sr^/a^       "^^12/^0  "  ^^2/^0 


+  Br,  e  +  Cr^  e 


] 


The  presence  of  excited  states  of  positronium  would  spread  out  the  wave 
function  and  one  would  therefore  expect  that  a  value  of  t  slightly  less 
than  -s  would  be  the  best  choice.  A  number  of  sets  of  values  of  the 
exponential  parameters  2,  s_,  t,  and  v  xrere  tried,  but  nothing  even 
approaching  a  systematic  determination  of  the  best  set  was  attempted. 
"Best*  is  here  defined  as  the  set  which  gives  the  lowest  upper  bound  on 
A.  Since  the  results  were  found  to  be  quite  sensitive  to  the  choice  of 
the  exponential  psrameters  the  realization  that  TT  satisfies  a  minimum 
principle  was  a  great  computational  advantage.  The  linear  parameters  were 
of  course  determined  variationally  for  each  choice  of  the  exponential 
parameters.  Results  are  given  in  Table  I. 

The  effect  of  polarization  is  rather  dramatic.   The  value  of  IT 
is  negative  so  that,  contrary  to  the  situation  in  the  static  approximation, 
the  positron  is  on  the  whole  attracted  to  the  hydrogen  atom.  Further, 
an  upper  bound  on  A"  now  yields  a  lower  bound  on  o'(k=o). 
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Table  I. 

Results  for  k  =  o,  with  the  trial  function  given  by  Eqs.    ()|.2l)  and 
(li.23)i   corresponding  to  the  set  of  exponential  parameters 
|q,s,t,vj  =  |o.^,  0.95,  0.U6,  O.78J  .     The  consistency  ratio,  R,  , 
is  defined  as  R.    ■  (A-  I  uilu  dt)  /A.     The  second  calculation  yielded 
a  trial  function  containing  no  nodes.     Note  that  while  the  second 
result  is  an  inprovement  over  the  first,   its  consistency  ratio 
deviates  further  from  xinity. 


A/a         B  CD       ;;pp^j       J°"^J 

o  bound  bound 

on  T/a  on  (^(k'o) 


\ 


-1.371ii        2.6516      -1.U620        0  -1.356  7.357na  ^      O.989 

0 

.l.Uli53        2.8106      -1.U979  -0.0532    -1.397  7.602na^^      0.966 
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V/e  note  incidentally  that  the  best  value  of  t  amon^  the  sets  considered 
was  t  =  0.U6,  i.e.,  somewhat  smaller  than  -^j  as  expected. 

ri2i 

Ilassey  and  Mcussa'-    -^   used  a  trial  function  of  the  fcnr. 
(ii.22)     rx{ryT^,r^^)  -  (2/a^^)^/^  e  "^^  *°  Tsin  kr^^ 


-r,/a  ~^i/*o 

+  (a  ♦  (b+cr^^g)  e  °)  (1  -  e    ^    °)  cos  kr. 


•]  ■ 


They  determined  "vT  variationally  for  various  values  of  k,   the  smallest 
of  which  was  k  =  0,2/a   .      Presiunably  they  avoided  k  =  o  because  they 
did  not  trust  the  fom  of  their  trial  function  for  that  value  of  k. 
If,   nevertheless,  we  use  their  form  we  find  X  ?5i  0,^12  a     ;  this 
value  does  not  differ  appreciably  from  that  obtained  under  the  static 
approximation.     It  should  be  noted   tliat  despite  the  extremely  poor 
value  for  A  the  abo'/e  calculation  provides  a  consistency  ratio 


K^   =  (A-U^^i 


dr)  /a  =  i.ou 


It  is  rather  surprising  that  the  two  trial  functions,  Eqs.  (h.21)  and 
(U.22),  each  containing  r, ^  dependence  and  three  variational  parameters, 
lead  to  such  strikingly  different  results,  corresponding  as  they  do 
to  an  effective  attraction  in  one  case  and  to  an  effective  repulsion  in 
the  other.  These  results  speak  strongly  for  the  need  for  obtaining 
bounds • 

It  may  well  be  that  the  possibility  of  the  virtual  formation  of 
positronium  makes  the  present  problem  particularly  sensitive  to  the 
choice  of  the  trial  function,  much  more  so  than  in  the  case  of  electron 
scattering  by  hydrogen,  for  example. 
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It  would  be  difficult  to  say  how  reliable  our  result  is  without 
obtaining  a   lower  bound  on  T.     We  have  been  unable  to  do   so,    for  we 
have  been  unsble  to  obtain  a   lower  bound  on  a    .      (The  difficulty  is 
connected  with  our  inability  to  prove  formally  that  a  composite 
bound  state  for  e"^  +  H  does  not  exist.)     We  therefore  sought  additional 
criteria  for  the  accuracy  of  our  results. 

In  particular,  since  no   composite  bound  state  exists,  we  exi:)ect 
by  analogy  with  the  one  body  problem  that  u  should  be  nodeless.     The 
trial  function  given  by  Eq.    (I4.2I)  has  nodes,   although  at  physically 
unint9r3sting  points    (r^  large).     Since  attempts  to   remedy  this  by 
varying  the  exponential  parameters  were  unsuccessful  another  term  \-ias 
added.     This   term  was  taken  to  be  proportional  to  exp(-5r^2/^o^'   ^^ 
one  thinks  of  this    term  as  being  related  to  the  virtual  formation  of 
positronium  in  its  first  excited  state  (it  does  not  have  precisely 
the  correct  form), one  might  expect  that  6  should  be  somewhere  near  lA. 
Further,   since  our  previous  exponential  factors  directly  coupled  r^ 
and  r^ ,  and  r„  and  r,pj  we  have  in  this  case  chosen  the  additional 
terra  to  be  proportional  to  exp(-2r,/a  ).     The  final  form  is 

(U.23)     u(r^,r2,r^2^  "  "  (Eq.(h.2l)) 

,,/  3a/2^     -^12/^^0 -V^o 

+  (2/a^-^)         Dr-j_e 

The  results  are  collected  in  Table  I. 

The  linear  parameters  were  determined  variationally,  i.e.,  without 
regard  to  the  question  of  nodes j  the  resultant  trial  function  is 
nevertheless  nodeless.  Further,  the  bound  on  J  was  lowered  only  slightly. 
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These  two  results  suggest  (in  the  absence  of  the  other  bound  on  A 
we  do  not  use  a  stronger  word)  that  our  bound  on  <f  may  well  be 

reasonably  clooe  to  the  true  value.  On  the  other  hand  our  trial  function 
does  not  indicate  the  presence  of  an  effective  polarization  potential 
which  falls  off  asymptotically  like  l/r  ,  as  demanded  by  perturbation 
theory '-■*, 

Our  trial  function,  Eq.  ([t.23),  may  now  be  used  to  estimate  the 
effective  range,  r  ,  for  the  e  H  system,  where  r  is  defined  by 


j'".^-^ 


(l;.2U)    r^  =  2  1  (u_^ -u^)  dr 


The  asymptotic  function,  u     ,  is  defined  by  Eq.   (3.7).     We  find 


(U.25)        r^  «i  Ua 
o  o 


This  very  large  value  of  the  effective  range  severely  restricts  the 
domain  of  applicability  of  the  shape  independent  approximation. 


(U.26)    k  cot^  «A  -I"""-  +  |k^r^ 


At  the  very  best  (if  the  coefficient  of  k^r      in  the  next  tern  is 
small*-    -"   and  if  the  expansion  is  in  fact  of  the  same  form  as  for  the 
one  body  case   ),  Eq.    (li.26)  should  be  useful  only  for  kr     somewhat 
less  than  unity,   i.e.,  ka     somewhat  less  than  0,1.     We  shall  return 
to  this  point  later. 


In  the  e  H  or  e~H  problem,  for  example,  "h,   does  not  vanish  with  k 
like  k      *      for  L50(due  to  the  long  range  polarization  potentialL-'J ). 
Here  the  effective  range  theory,  in  its  usual  form,  would  not  be  applicable, 
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It  might  well  have  been  less  time-consuming,   for  all  of  our 
calculstionsj   to  have  fixftd  the  exponential  parameters  arbitrarily 
and  to  have  used  more  terms   containing  linear  parameters.     The  present 
approach  was  used  in  orxier  tc  have  a  good  trial  function  in  reasonably 
simple  form  for  other  applications;   in  particular,  we  will  now 
consider  a  calculation  at  a  non-vanishing  energy  with  a  trial  function 
which  is  a  direct  extention  of  the  zero  energy  form  given  in  Eq.    (U.23), 

C.   Polarization  Considered,  k  j^  o 

In  order  to  shoxir  that  p^  =   00  without  introducing  cut-off 

potentials  we  have  restricted  ourselves  to  the  special  case  of  zero 

energy  scattering.    It  is  not  to  be  thought,  however,   that  calculations 

leading  to  a  bound  on  the  cross  section  can  be  performed  only  at  zero 

energy.     As  indicated  in  Section  3>  for  sufficiently  small  non-zero 

energies  the  cut-off  point  may  be  made  large  enough  so  that  the 

difference  between  rj  and  vj"     (the  phase  shift  in  the  presence  of  the 

cut-off  potentials)   is  negligibly  small.     To  illustrate  the  procedure 

we  describe  here  a  determination  of  a  loirer  bound  on  Y)     for  ka    =    0.2. 

(c 

To  begin  with  we  ignore  the  question  of  obtaining  a  bound  and 

perform  an  ordinary  variational  calculation  for  the  true  problem.  We 

use  a  trial  function  which  reduces,  as  k  ->  0,  to  the  zero  energy 

function,  Eq.  (i|.23),  namely 


(U 


.27)   ^/2V(^1'^2'^12^'^  "   (2Ao-^)^'^^[Ae"'"^  ^°  (l-e"'^  ^°)  cos  kr^ 

-  e      sm  kr,  A  +  Br-,e 


-tr^2/%  -  ^^2/^0   ^   -^12/^^0  -  2^1/^0 


+  Gr-e  +  Dr-|_e 


]  ■ 
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From  the  asymptotic  fom  of  this  function  we  see  that  A  =  -tanlTA- 
The  variational  principle  is  now 

r 

(U.28)     k  tan-h  ^  k  tan>i  +       ^/2\n  "^  "6./^     ^'^ 

where cC  is   given  by  Eq.    (2.2).     The  result  of  the  calculation  is 
T]      =    0.l6a5     , 

^     =J  ^4^^/2>^^    =     -0.0017/a^ 
and  11  5^  0.156.     As  a  point  of  comparison  we  note  that  a  variational 
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calculation  (including  polarization)  performed  by  Massey  and  Moussa'-    -" 
for  the  same  scattering  energy  gave  77  7^  -  O.O98  • 

Since  the  trial  function  is  based  on  the  zero  energy  function  which 
we  expect  to  be  quite  good,   we  have  obtained  what  we  believe  to  be  an 
accurate  estimate  of  TJ  ,  although  not  necessarily  a  lower  bound. 
Before  considering  the  calculation  involving  the   cut-off  potentials 
we  recall  that  the  cut-off  point,  a  ,  is  limited  by  the  condition  ka  <  n-9, 
In  order  to  be  able  to  pick  the  largest  value  of  a  (for  fixed  k)  we 
wish  to  choose  the  smallest  value  of  6  consistent  with  the  condition 
TJ  <  &.      (Note  that  we  desire  that  this  relation  be  satisfied,   in  order 
that  our  bound  be  a  useful  one.     As  discussed  previously,   the  bound  will 
be  valid  in  either  case.)     The  trial  function,  Eq.    (l4.27)>   is  noimalized 
with  ©  "  (1/2V1.     This  would  lead  to  a  ■  7.9a  .     In  order  to  obtain  a 
variational  estimate  of  f)  corresponding  to  the  use  of  a  trial  function 
with  different  (and  smaller)  9  -  nonnalization,  we  may  simply  choose 


(a. 29) 


^  =  ^/2y  °osysin(|-9) 
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so  that 

(U.30)         k  cot   ffi-e)A::k  cot(7j-9)  -     cos  19 /sin ("n  -  9) 


I 


where  Tfl  and  I  have  already  been  determined.     Vie  expect  that  this 
variational  estimate  of  TJ"  should  change  very  slowly  as  9  is  decreased 
from(l/2)n.     It  is  found,   in  fact,    that  to  three  significant  figures 
the  result  is  unchanged  for  9  as  small  as  O.3U. (Smaller  values  of  9 
lead  to  appreciably  lower  estLmates  of  tT.  )     We  have  therefore  chosen 
9  =  0.3U00  which  allows  us   to  pick  a  =  liia   • 

To  obtain  a  bound  on"^     we  must  now  perform  a  similar  calculation 
with  oo  replaced  by  Aj     which  contains  the  cut-off  potentials,  Eq.    (3.5)» 
The  trial  function,  vu      ,   must  now  satisfy  the  additional  requirement 
that  Jo  Ug       =  0   for  r-,   >  a",  while  lu       must  of  course  be  continuous 
in  slope  ana  value.     Such  a  function  may  be  written  as 

('4.31)     u^^pfcos  -y/sinCy-d)!       =  iUy2V  (Eq.   (U.27))     g(r^) 


+   (2/a^3)l/2^_g(^_^))  ^""2  \in(kr3^+|)/cosTj 


where 


(a.32) 


g(o)  =  1 
g'(a)  =  o 
g(r^)  =0  ,  r^  >  a 


If  g(r, )  is  close  to  unity  throughout  the  region  r,  <  a  we  might 

expect  that  the  results  obtained  will  not  differ  very  much  from  those  of 


This  may  be  seen  from  the  fact  that  In  order  to  obtain  Pa  =  00,  p  must 
v.anish  beyond  r^  =  a.  However,  J  (JCt.  )^p'''^dV    should  exist. 
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the  above  variational  calculation.  With  this  in  mind  we  chose 


(U.33)  g^^i)  = 


1  -  e 


-  ^(1-r^/a) 


1  -  e 


< 


r,  <  a 


r,  >  a 


where  ^  is  arbitrary.  If  we  let  Y,   take  on  a  Isr^e,  positive  value 
fe  «i' 5a,  say),  g(r^)  will  have  the  desired  shape. 

We  may  write  ,  ' . . 


(U.3U)    "a  .<^o  "a  .  ^'i 
w,c  c  w,c 


y 


u^^u^dl?  +7^  . 


Due  to  the  shape  of  g(r-),  contributions  to  the  Integra]  "J?  are  very- 
small  for  r^  <  a  }  for  Pj^^  a  the  expression  is  dominated  by  exponential 
factors  which,  owing  to  the  large  value. of  a,  are  also  very  small.  It 
is  in  fact  easily  determined  that  to  the  accuracy  of  our  calculation"!^ 
is  small  enough  to  be  neglected  entirely.  Therefore,  the  variational 
estimate  of  to  is  actually  a  bound  on  ^  and  hence  is  extremely  likely 
to  be  a  bound  on  Ti  itself,  since  the  contribution  of  the  potential  in 
the  region  r-  >  a  is  expected  to  be  extremely  small.  Results  appear 
in  Table  II. 

It  should  be  mentioned  that  even  neglecting  the  difference  between 
Ti  and  "Jj  we  still  do  not,  contrary  to  the  zero  energy  case,  obtain  a 
rigorous  bound  on  the  cross  section,  defined  by 


(a. 35)    cr(k)  =  (hnA  )  sin 


The  substitution  of  y^ower  *  ^'^®  calculated  lower  bound  on  the  phase 
shift,  into  Eq.  (li.35)  will  yield  a  lower  bound  on  cl'(k)  provided 
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Table  II 

Results  for  ka     =  0.2,  with  the  trial  function  given  by  Eqs.  (U.27), 
(I1.3I)  and  (U.33).     The  exponential  parameters  in  Eq,  (U*27)  are  the 
same  as  those  cited  in  Table  I.     The  bo\ind  is  obtained  onyj  ,  the 
phase  shift  for  the  problem  in  which  the  positron  interaction  vanishes 

outside  a  sphere  of  radius  lUa   .     y\  variational  calculation  by  Massey 

Till  — 

and  Moussa,  *■    -^ ,   for  the  same  energy,   gave -nfd  -  O.O98. 


A /a  B  C  D  Lower  bound 


0„| 


c 


-  0.8285     2.0U3U     -  1.1271       0  0.1$b 

-  0.8300      2.1159      -  1.1U6U      -0.0lil2  0.156 
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-n  <  n-yj^  .      Such  an  upper  bound  has  not  been  obtained;  however j 

in  vieTv'  of  the  expected  accuracy  of  our  variational   calculation,    it  is 
almost  surely  valid.     The  difficulty  does  not  arise   in  the  limiting 
case  of  zero  energy  since  there  the  corresponding  piece  of  information 
is  supplied  by  the  knowledge  that  no   composite  e  II  bound  state  exists. 

It  is  easily  shown,   in  the  one  body  problem,   that   if  the  potential 
vanishes  beyond  r  -  a,   and  further,   if  ka  +  "rv  <  n  ,   then  the  wave 
function  is  nodeless  in  the  region  r  <  a.      Thus,   the  standar-d  definition 
of  the   phase  shift,    i.e., 


(I1.36)        •>)   =  lim      (niT  -  kp, 
*      n->oo  \ 


can  be  written,  in  this  casa,  as 

(U.36')  t[=  nn  -  kp^   ,         ?„  >  »    » 

where  p     is  the  nth  zero   (not  including  the   zero  at  the  origin)   in  the 
wave  function.      Now  consider  the  point  r'    =   (ii-'h)A:.     From  the 
assumed  restriction  onil   it   is   clear  that  r     >  a.     Now  the  wave 
function  takes  en  its  asymptotic  fonn  beyond  the  point  _a,   and  must 
therefore  vanish  at  r'.     But  from  Eq.    (U.36),   r'   =  p,    ,   so  that  we 
have  shown  that  the  fii^^t  zero-  occurs  beyond  a.     In  the  present  problem, 
under  similar  conditions,    (i.e. ,  W.   +  W_^    =0   for  r,   >  a,   and 
ka  +  T7  <  n)   it  may  be  shown  that  the  function 


\ 


n 


2 


(li.37)       g(r-j^)  5         T^     dr2  dp  V.{r^)  u(r^,r2,p) 


J 

is  nodeless  for  r,  <  a.  For  s_  large  enough  it  seem.s  reasonable  (we 
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have  not  proved  it)  that  a  stronger  ststement  can  be  made,  namely, 
that  u  itself  is  nodeless  for  r^  <  a.  If  we  assume  this  statement 
to  be  true,  it  may  be  used  to  test  our  trial  function.  With  D  =  o 
the  function  has  nodes  for  r,  <  a  but  these  nodes  do  in  fact  disappear 
when  the  full  trial  function  is  used.  Further,  this  im.provem.ent  in 
the  trial  function  is  accompanied  by  only  a  small  improvement  in  the 
bound  on  the  phase  shift.  Both  of  these  considerations  suggest  that 
this  bound  may  well  be  in  the  neighborhood  of  the  true  value. 

We  have  previously  noted  that  the  shape  independent  approximation 
can  at  best  be  expected  to  be  applicable  for  ka  <  0,1.  If, 
nevertheless,  we  apply  Eq.  ()i.26)  to  the  present  case,  ka^  =  0.2, 
we  find  tan  ^  f^  0.21.  This  differs  considerably  from  the 
variational  result,  tan  ^-^0,16.  The  difference  is  ven--  likely  due 
to  the  inapplicability  of  Eq.  (li.26),  but  we  may  not  exclude  the 
possibility  that  our  calculated  value  of  r  is  too  large  due  to 
inaccuracies  in  our  zero  energj'-  wave  function.  It  is  of  course 
also  possible,  though  less  likely,  that  our  variationally  calculated 
values  of  tan  ^A  for  k  =  o  and  0.2  are  not  close  enough  to  the 
correct  values. 
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6.  Discussion 

We  have  shown,  by  providing  a  particular  example,  that  it  is 

not  only  possible  to  obtain  a  bound  on  the  cross  section  for  scattering 

hj   a  compound  system  but,  under  certain  circumstances,  that  no  additional 

work  beyond  that  of  performing  a  variational  calculation  is  required. 

This  is  true  in  the  case  of  zero  energy  scattering  when  no  composite 

bound  state  exists.  For  slightly  larger  energies  the  bound  obtained 

is  no  longer  rigorous,  since  the  interaction  between  the  scattered 

particle  and  the  scattering  system  is  neglected  outside  a  certain 

region,  VJe  believe,  however,  that  in  the  example  presented  above, 

where  the  cut-off  point  is  at  llta  ,  the  neglected  potential  contributes 

practically  nothing  to  the  phase  shift  and  that  our  bound  is  extremely 

likely  to  be  valid. 

The  fact  that  the  variational  principle  for  cot  Tj"  may,  in  special 
cases,  be  a  minimum  principle  as  well  has  some  interesting  consequences. 
Generalizing  the  results  obtained  previously,  we  may  state  that  the 
scattering  length  for  positron  scattering  from  a  static  atom  is  greater 
than  the  true  scattering  length.  A,  provided  that  the  positron  and  the 
atom  do  not  form  a  bound  state.  Also,  a  variational  estimate  of  A 
will  necessarily  be  positive  in  this  case  unless  the  trial  function 
contains  some  dependence  on  the  separation  between  tlie  positron  and  one 
or  more  of  the  electrons.  (This  is  analogous  to  a  result  obtained  by 
Ore'-  -^  in  his  investigation  of  the  stability  of  simple  positron  compounds.) 
Another  result  of  this  formalism  which  would  be  otherwise  rather  difficult 
to  prove  is  that  under  the  conditions  for  which  p/.  /A  =  oo  the  Kohn 
form  of  the  variational  principle  for  the  phase  shift  (derivable  from  the 
Kato  form  by  setting  9  =  n/2)  has  at  least  one  important  advantage  over 
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the  Hulthen  formulation.  That  is,  an  increase  in  the  number  of 
parameters  in  the  trial  function  guarantees  an  improved  result. 
This  remark  applies,  for  example,  to  the  problem  of  zero  energy  triplet 
scattering  of  electrons  from  hydrogen. 

The  simple  result  that  at  zero  energy  the  ordinary  variational 
calculation  provides  a  bound  if  the  system  has  no  composite  bound 
state  suggests  that  the  Kato  mechanism  is  not  necessary  in  this  case, 
since  the  above  conditions  make  no  reference  to  the  associated  eigenvalue 
problem.  This  is  easily  verified  by  observing  that  in  the  absence  of 
a  bound  state  the  operatorcC  is  negative  definite  in  the  space  of 
quadricatically  integrable  functions.   Now  for  0  <  &  <  n. 


(5.1)     w^  =  lim  w^A 
k*o 

has  the  form  of  a  zero  energy  bound  state  wove  function.     That  is 


(5.2)  iL  ->  const.   R  (rj)  »  for  r-j^  ->  oo. 


Thus  w~  may  be  considered  to  be  the  limit  of  a  sequence  of  functions. 


Wo(e),  which  have  the  asymptotic  form 


e 


-cr, 
(5.3)     ^Q^^)  ">  const.  ^  (r2)  e     ,      for  r^  ->  oo  . 


Since 

(5.U) 


r 


Wa(e)Xwa(e)dr     <    0 


"■9""'^  ™e 


we  conclude  that  this  is  true  for  e  =  o  as  well.  This,  along  with 
Eq.  (2.6),  gives  us  the  desired  result.  Note  that 

Wq  ->  const.  R(r2)r^  ,  for  large  r^, 
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so  that  the  theorem  is  not  true  in  general  for  this  choice  of  9  -  normalization. 
We  had  also  come  to  this  conclusion  from  the  point  of  view  of  the  Kato 
formalism.  (See  Section  3). 

While  this  proof  merely  reproduces  a  result  previously  obtained 
it  is  still  quite  interesting.  First,  our  conclusions  are  freed  from  ' 
any  dependence  on  the  validity  of  the  generalization  of  Levinson's 
theorem  (see  Section  3).  More  significantly,  the  type  of  reasoning 
involved  points  directly  to  the  extension  which  includes  situations 
in  which  bound  states  do  exist.  Here  the  results  we  have  obtained  do 
net  follow,  in  any  simple  way  at  least,  from  the  Kato  formalism  and 
may  provide  a  considerable  improvem.ent .  Furthermore,  we  have  been  led 
to  the  development  of  a  method,  valid  at  higher  energies,  which  reproduces 
all  results  of  Kato  for  which  p_=  oo  and  which  should  otherwise  provide 
a  useful  supplement  to  the  Kato  method.  This  work  will  be  reported  on 
in  the  future. 
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Appendix  A 


In  the  numerical  evaluatio 
integrals  of  the  form 
(A.l) 


n  of   I  ujb  u 


dti  (Section  UB)  various 


00  /OO  /'I'T, 

l(a,p,Yj  i,m,n)  2  ]       dr^   1       drj 


7 


A^^2  ^r2-pr^-rr^2     Z      m       n 

dr^2^  ^2  ^1  ^12 


l^rgl 


occur.  These  integrals  may  be  generated  from  the  basic  integral 


(A.2)     I{a,?,r;   0,0,0)  «  (a^^)(0!yJ(a4Y)  =  | 


using  symmetry  relations  such  as 


(A. 3)     I(a,p,Yj;2,m,n)  =  I(p,a,Y}  m,J?,n) 


and  recursion  relations  such  as 


(A.h)  l(a,p,Y;  JP+l,m,n)  =  -  -^  I(«i,P,Yj/,m,n) 


and 


r 


(A. 5)  l(a,p,Yj/-l,m,n)  =  I(a',|i,Yj  i',m,n)  da' 


00 


/ 


We  list  below  expressions  for  those  integrals,   involving  non-vanishing 
values  of  y,  which  entered  into  the  calculation. 
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(A.6) 


l(a,n,Y;  1,0,0)=  I  [^-^] 

l(a,%Y;   1,1,0)  =  -^       ~  +  a  +  p   +  Y 

1 

I(a,-,Y;   0,1,1)   =  li»       -1.  +  a   +  p,   +  Y  J 


l(a,3,Y;    1,0,1)   =li^|-|-  +  a   +  B+Y 


7   l"^ 


i(a,.,Y;i,i,i)=^[2.^.:^.^] 
i(a,?,Y;  1,-1,1)=  -  -  j-^r^  lo 


(^2:^  --  p^   T^V? 


In  the  calculation  involving  the  non-vanisbj-n^  scattering  energy 
(Section  hC)   additional  integrals  had  to  be  evaluated.      The  results 
nay  be  stated  in  terms  of  the  general  form 

(A. 7)     l(a,p,Y;  i,m,n;   p,q)  = 


i 


00       ^00       (^i^^2 


dr. 


dr. 


^ 


dr^2  e-^^2-Pri-Yri2  rf  r^^^g"  cosPkr^sinV^     . 


VrgI 


Again  we  quote  only  those  of  the  required  integrals  for  which  y  j^  0, 


(A. 8)  l(a,?,Y;  1,0,0;  0,1)  =  2k\       ^^%^  | i^— ^      ^ 


(a^-r^)^  L(P+Y)^+k^   (a+p)W 


, 1_  r    2(a-^p)    "I 

'  a^-Y^  L((a+p)W)^  J 
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I(a,p,Y;  1,-1,1;  0,1)  =  iaWr'l^  tan^  ^  -J^^ 


+  a  <- 


I(c,S,y;  1,0,1,1,0)  =4  l(a,P,YJ  1,-1,1;  0,1) 


ok 


/   2     2n-2     l6aY  a+B 

(a   -Y  )      <    b   U   •  p — ? 

a^^-^C         (a+p)^+k 


Ua  8c  k^ 


(P+Y)^+k^        [(P  +  f)Wj 


+     a  <- 


I(a,p,Y;   1,0,0;  1,0)  =  -^   Tt^Vt  "  r^rH  1 

(a^-Y'^)      L  (P+y)   +k^        (a+p)^+k^J 


(a^-Y'^)      L  (P+y)   +k^        (^+P 


2        I         1  2(a+l3)^ 


a  -y'    L(a+p)'^+k^        [fn.n^^.v.^V 


r(a+p)Wl' 


8aY  ,      //„.    ^2.,2. 


l(a,?,Y;  1,-1,1;  1,0)  =       J^T  o   •  log((G+Y)  +k  ) 

,        2a  ^        2(P+y) 


+  a  <- 


l(a,.P,Y;   1,0,1;  0,1)   =  -  A  l(a,p,Y;   1,-1,1;   1,0) 


«     l6kaY       ,  __1 +       8ka  ^         pn* 


2  ,  2  ■*■  772     2T? 


(aV)^       (P^)W       (aV)^       [(p^)2,,2-j 


+  a  < >  Y 
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Apperdix  B 

We  wish  to  consider  here  a  point  which  srises  in  the  generalization 
of  the  Kato  formalism  from  the  case  of  scattering  by  a  center  of  force 
to  the  situation  in  which  the  scatterer  is  a  compound  system.  Suppose 
we  are  dealing  with  a  problem  of  the  type  discxissed  in  Section  2. 
Then,  according  to  the  definition  of  the  eigenfunctions  9'  g  j 
(solutions  of  Eq.  (2.9)  with  boundary  conditions  given  by  Eqs.  (2,11)), 
one  can  easily  construct  a  scattering  problem  for  which  the  set  of 
functions,  Y  0  ^  I ,  is  not  complete.  Thus,  with  W.  =  0  and  p  =  pir^) , 
the  eigenfunctions  will  each  be  of  the  form 

The  tine  solution  to  the  real  problem,  u„  ,  will  also  be  of  the  product 
form.  If  one  chooses  a  trial  function,  Ug.   ,  which  is  not  simply  a 
product  of  R(r  )  with  some  function  of  r,  then  w  =  u^  -  u-  will  not 
be  of  the  product  form  either  and  clearly  can  not  be  represented  as  an 
expansion  in  the  9^  q»  Of  course,  this  particular  case  represents  no 
real  difficulty  since  one  can  easily  recognize  the  above-mentioned 
restriction  on  the  form  of  9^  „   and  u_  and  choose  a  trial  function  of 
the  same  form.  It  is  possible,  however,  that  the  true  solutions  may 
satisfy  some  symmetry  properties  which  are  not  easily  recognizable. 
It  would  be  desirable,  therefore,  to  have  a  precise  formulation  of  the 
conditions  under  which  the  set  j  ^j.  q  i  is  complete.  We  shall  not 
attempt  to  provide  this  here.  However,  some  general  conjectures  along 
these  lines  may  be  made. 
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We  obserre  first  that  the  boimdary  conditions  placed  on  the 
eigenfunctions  in  the  one  body  problem,  i«e. , 


(B.l)  0     ^(r)  — *-     const.   sin(kr  +  9  +  nn)   ,       as  r, — >■     oo 


0^^^(C)       -     0 


are  ccmpletely  equivalent  to  the  requirement  that 

'00  y-OO 

(B.2) 


/ 


^tJ"  ^^  o  <^r  -         (^„  Q  oC  Wq  dr  =  0 
w      TifV  J      n,w        y 


wheiB  «L  is  now  the  appropriate  operator  for  the  one  body  problem. 
Thus,  a  solution  of  the  equation 

<£^  +  HPfZf  =  0 
has  the  as3nnptotic  form 

0  — >     C  sin  (kr  +  0  )  ,         as  r-^  <x> 

where  C  is  a  non-vanishing  constant  (i«e,,   there  ar«  no  bound  states 
in  the  continuum).     We  may  therefore  write 

rOO  /<X)  00 

/\  3    j       w^<C  0  dr  -  J     0<Cw^  dr  -  (w^i/'  -   0w^)   | 
o  'b 

»  j6(0)  w^(0)  +  kCD  sin(9  -(£>) 
using  the  fact  that  w_(0)   «  0  and  w„— *•  D  sin(kr+©)  as  r --»•  oo,  D 

I 

being  a  non-vanishing  constant.     We  note  that  w„(0)  is  an  arbitrary 
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number  depending  on  the  particular  choice  of  the  trial  function. 
Since  CD  /  0  we  must  therefore  have  0(0)  =0  and  @  =  C^  in  order 
that  j/^  vanish  for  an  arbitrary  function  Vr.  satisfj'lng  the  boundary 
conditions.     The  equivalence  of  Eqs.    (B.l)  and   (B.2)  is  then  apparent. 
The   generalizations  of  these  equations  in  the  many  body  problem  are, 
however,   not  necessarily  equivalent.     To  see  this  we  need  only  return 
to  the  problem  cited  above,   in  which  W       =  0.     For  certain  values 
of  ij.  =  u     there  exist  solutions  of  the  form 

where  R     j^r^)  represents  an  excited  state  for  particle  2  and  F     ,(r, ) 
ni,L     2  ^  m,L     1 

is   a  bound  state  function  (vanishing  asymptotically).     We  note  that  this 
type  of  solution  will  satisfy  the  condition 


j   w/0^c 


r 


(B.3)  I   w^<)o0_dt  - 


'^rfi%  dt 


since  the  surface  terms  vanish j  however,  they  do  not  satisfy  Eqs.  (2.11). 

- 

Further,  the  inclusion  of  these  solutions  in  the  set  -  0  ,,  •  allows  the 

n,w 

possibility  that  the  set  will  be  complete  even  for  trial  functions  not 
of  the  product  form.  It  would  therefore  seem  to  be  preferable,  in 
general,  to  define  the  eigenfunctions  0     as  those  solutions  of  Eq.  (2.9) 
which  satisfy  Eq.  (B.3).  The  question  of  whether  this  is  sufficient 
to  ensure  completeness  would  still  have  to  be  studied,  of  course. 
Further,  even  assuming  completeness,  one  would  then  have  to  investigate 
what  effect,  if  any,  this  more  general  definition  of  the  set  1  )^  V  had 
on  the  determination  of  the  eigenvalues  a^  and  p^.  These  considerations. 
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however,  do  not  alter  the  validity  of  our  present  results  for  the 
e*H  problem  and  it  is  hoped  that  in  general  they  introduce  no  major 
obstacles  to  the  practical  application  of  the  Kato  method.     This  would 
be  the  case  if,   for  example,   the  non-decaying  eigenfunctions  were 
complete  for  a  well  defined  class  of  functions  to  which  w^  belongs. 
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